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Abstract 



We consider the spectral form factor of random unitary matrices as well as 

of Floquet matrices of kicked tops, as given by the (squared moduli of) the 

traces t„ = TrF" with the integer "time" n = 0, ibl,ib2, .... For a typical 

,^ , matrix F the time dependence of the form factor |tnP looks erratic; only after 

(N| I a local time average over a suitably large time window An does a systematic 

^^ I time dependence become manifest. For matrices drawn from the circular 

(•^ ■ unitary ensemble we prove ergodicity: In the limits of large matrix dimension 

^^ • and time window An the local time average has vanishingly small ensemble 

fluctuations and may be identified with the ensemble average. By numerically 

diagonalizing Floquet matrices of kicked tops with a globally chaotic classical 

'^ ' limit we find the same ergodicity. As a byproduct we find that the traces 

Q ■ in of random matrices from the circular ensembles behave very much like 

2 ' independent Gaussian random numbers. Again, Floquet matrices of chaotic 

Q . tops share that universal behavior. It becomes clear that the form factor of 

chaotic dynamical systems can be fully faithful to random-matrix theory, not 

is^ ^ only in its locally time-averaged systematic time dependence but also in its 

Vh ' fluctuations. 
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I. INTRODUCTION 

In the present paper we propose to discuss the time dependence of the form factor of 
individual dynamical systems and of individual random matrices and compare with the 
ensemble based predictions of random-matrix theory. Our investigation was motivated by 
recent efforts to extract universality from semiclassical arguments. Among these, Smilansky 
and coworkers |^J^ have argued that random-matrix behavior should prevail in the autocor- 
relation function of the spectral determinant of dynamical systems, at least after a suitable 
average; the conclusions were based on assuming Gaussian statistics for the form factor 
fluctuations. Similarly motivated work by Prange and Kunz |Q] has partial overlap with 
ours. 

In order to look for universal spectral fluctuations from a given quasi-energy spectrum 
one must first "unfold" the spectrum to constant mean level spacing, or else systematic 
system specific variations of the density of levels across the spectrum would be mixed up 
with potentially universal fluctuations. Moreover, smoothing over certain evergy intervals, 
denoted by the overbar in p(E + AE) p{E), for the correlation functions of the density of 
levels or over suitable time intervals for its Fourier transform, the form factor, is necessary. 
Without such averaging the form factor would display erratic variations not at all resembling 
the smooth behavior predicted by random-matrix theory for averages of matrix ensembles. 
The level density correlator would even consist of delta peaks. 

If a conservative quantum system with global chaos in its classical limit displays universal 
spectral fluctuations a la random-matrix theory, such universality may reveal itself in the 
distribution P{s) of nearest-neighbor spacings, for which random-matrix theory predicts the 
power law P{s) ~ s'^ as s — >■ 0, and the degree P of level repulsion is determined by symmetry, 
most importantly the absence or presence of time reversal invariance. Equally popular and 
even more easily amenable to theoretical reasoning is the two-point function of the level 
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density C{AE) = p(E + AE) p{E) — p{E) , for which universal behavior implies the same 
power law for small energy difference AE. The Fourier transform of that correlation function, 
the so-called form factor K{t), may of course also be used as an indicator of universality 
(or of lack thereof). Ergodicity of the underlying classical dynamics |^ imparts a linear 
dependence on the time r, K{t) ~ |r|, for |r| larger than the shortest periods of classical 
periodic orbits, but smaller than the Heisenberg time th (the time scale corresponding to 
the mean level spacing as a unit of energy). For larger times, \r\ 3> th, the form factor 
tends to a constant value, and the character of the approach is related to the degree of level 
repulsion. 

For the sake of concreteness we base our work on unitary N x N matrices, such as 
arise as Floquet matrices F for periodically driven systems with a compact phase space and 
thus finite dimensional Hilbert space. As examples of such dynamics we shall take kicked 
tops, ones with the spherical phase space dominated by chaos and, on the other hand, 
symmetries chosen so as to fit any of the three principal universality classes (orthogonal, 
unitary and symplectic). The corresponding random-matrix ensembles are, of course, the 
so-called circular ones of Wigner and Dyson (COE, CUE, CSE). 



Our findings within random-matrix theory extend the ones previously presented in |TT 
and may be summarized as follows. They all concern the "traces" t„ = trF" with n = 
1,2,..., which give the form factor as K{n) = |tnp; the integer exponent n serves as a 



discrete time, counting the number of periods of the external driving. (In the large-A^ 
hmit we shall eventually introduce r = n/N as a quasi-continuous time.) The circular 
ensembles yield marginal probability densities {5{t — tn)) which, in the limit of large matrix 
dimensions, A^ — > oo, assign Gaussian statistics to finite-order moments. In particular, 
the first four moments bear relations of Gaussian character, as if coming from P{t) = 
(7r|tnp)~^ exp f — |tp/|t„|2j. For different traces t„ and tm we show, again for the large- 
A^ limit, that the unitary circular ensemble does not give cross-correlations, in the sense 
l^npl^mp/d^nP |^mp) — 1 = 0{1/N) ioT 771 ^ u. Usiug thcsc rcsults for ensemble averages 
we show for the CUE that the form factor is ergodic in the large- A^ limit: The time average 
(l^nP/l^nP) = (^^)~^Sn't^" l^n'P/l^nP ^^^ an eusemblc variance vanishing in the limit of 
a large temporal window An as 1/An. This means that with overwhelming probability 
every random unitary matrix drawn from the appropriate circular ensemble has the same 
time averaged form factor and that the latter equals the ensemble-averaged form factor. 
Inasmuch as a dynamical system has a Floquet matrix typical for the appropriate ensemble 
one can expect universality for its time averaged form factor as well. 

By numerically diagonalizing Floquet matrices of kicked tops from the various universal- 
ity classes we have calculated, for each of these, the form factor K{n) = |t„p and its time 
average over a finite window An. Normalizing to the ensemble- averaged form factor |t„p a 
la random- matrix theory we throw all |t„p/|trip for a given top into a histogram and find 
this to reproduce the Gaussian behavior predicted by random-matrix theory. The ergodic 
character of random matrices is also respected in full by the Floquet matrices of chaotic 
tops: The time-averages (|t„P/|t„|2) come out to have variances within the respective data 
sets {tn} varying with the time window An as 1/An. 

II. THE FORM FACTOR AND ITS FLUCTUATIONS IN ENSEMBLES OF 

RANDOM MATRICES 

The density of eigenvalues e~"^* of a unitary matrix F can be written as 

-\ N oo -| oo 

^^^' 2nN ^, ^ 2-kN ^ ^ ' 

%=\n= — OQ i"ii> n= — OQ 

We herein meet the traces 

AT 

^n = Ee"'"'''=tr^", n = 0,±l,±2,... . (2) 

i=\ 

If F is drawn from a homogeneous ensemble of random-matrices like any of the familiar 
circular ones of Wigner and Dyson, one has the ensemble average tntn' = ^n,-n'|^np and 
thus the two-point correlation function 



C(e) = p(^_+e27r/Ar)p(^) _ f (3) 

_ V^ U |2„i27rne/7V 
/ y |''n| e 

The mean form factor |t„p is well known for the circular ensembles |I^pT]] and will be given 



further below. Since we want to determine higher-order moments like |tn|^; it is best to 
consider the characteristic function 



Pl^j^{k) = exp — ^(A;e-i"'^> + k*e''"^*) (4) 



N 



Y[exp{-i\k\ cos{n(pi)) 



i=l 



where the joint densities of eigenphases of the circular ensembles [0,|ril] may be used to 
calculate the averages for (3=1 (Circular Orthogonal Ensemble), (3 = 2 (Circular Unitary 
Ensemble) and (3 = 4 (Circular Symplectic Ensemble). The homogeneity of these ensembles 
entails the characteristic function to depend on k and k* only through the modulus \k\ which 
we shall simply denote by k henceforth. It follows that all odd-order moments vanish, as 
well as those even-order ones where a trace t„ is not accompanied by its complex conjugate 
t* to form powers of |triP- The non- vanishing moments are 

\^=(-^rM '^pUk)\.^o. (5) 



III. CIRCULAR UNITARY ENSEMBLE 

The characteristic function is easily evaluated using any of the standard methods of 
random- matrix theory and takes the form of a Toeplitz determinant, 

oo 

P^j,{k) = detM, M^,= ^ (-i)l^lj|,|(A;)5^,,_„, , (6) 

s=— oo 

with Js the Bessel function of integer order. The desired Taylor coefficients of det M{k) are 
most conveniently evaluated through the expansions 

M{k) = 1 + A;M(i) + — M(2) + :^M(3) + . . . , (7) 

det M{k) = e'^ '-''^^ = exp(^r2 + ^t, + ...) . 

The derivatives M*^*) = d*M(A;)/dA;*|fc=o can be calculated with the help of Js(0) = 5^,0 and 
dJ,(fc)/dA; = (J._i(A;) - J.+i(A;))/2 as 

<.) = -^(<5;-" + <5;+"), (8) 



i S? 
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Inasmuch as the Taylor coefficients of det M(fc) are related to the moments Itn^™ and the 
Taylor coefficients of its logarithm In det M = tr In M to the corresponding cumulants, their 
interrelations 



n = tr M(^) , (9) 

r2=tr(M(2)-M«2), 

rg = tr (M(3) + 2M(^)3 _ 3^(1) M^^)) , 

r4 = tr (M(^) + 12M(^)2^(2) _ q^(i)4 _ ^m^2)2 _ ^m^D m^^)) 

resemble the familiar ones between moments and cumulants of random variables. One arrives 
at Ti = rs = and 

{ -n/2 Q<n<N 

^' = \-Nl2 N<n , (10) 

( 0<n< N/2 

r4 = i |iV - |n N/2<n<N (11) 

[ -|A^ N<n 

and these give the well-known CUE form factor 



n < n < N 



l*^"' ~] N N<n ^^^'' 



as well as its variance 



Var(|t„n = |t„|4 - |t„|2 =4r^ + -r^ (13) 

n^ 0<n< N/2 

n^ -2n + N N/2 <n<N (14) 

N^-N N <n 



|tnP +\t2n\^-2\tn\^. (15) 



Within the interval < n < N/2 the mean |t„P and the variance Var(|t„p) are related as if 
tn had a Gaussian distribution 

P(t„) = ^e-l*"l'/^ , (16) 

" I ■'ml 

and that distribution also correctly reflects the vanishing of all odd-order moments and 



cumulants. In ||Tl[ we had established the more general result that ([T6| ) even reproduces all 



moments |t„|2™ with n < N/m. Below we shall show that the Gaussian relation between 
\tn\'^ and Var(|t„p) prevails even for n > N/2, up to asymptotically negligible corrections of 
order ln(A^)/A^. 



IV. CIRCULAR ORTHOGONAL AND SYMPLECTIC ENSEMBLE 

For the COE and CSE the evaluation of P^^(/c) in (§) proceeds analogously ll| and yields 
the Pfaffians PJ^Nik) = JdetAP{k), 



^k'(^) = E -Mik) J\s'\{k) ^SIZ^ <5^^-«(^+^') , (17) 

I lib \~ I tj 

s 

wherein the underhned indices m,rr]! are semi-integer and in the ranges \m\, |m'| < {N — l)/2 
for /3 = 1 and |2?i|, \rn/\ < {2N—l)/2 for (3 = 4. One encounters more comphcated expressions 
for the expansion of ln74^(A;) since A^{0) ^ 1 and an additional factor 1/2 in the ti due to 
the square root in the Pfaffian. To present our results in a concise form we introduce the 
shorthands 

^1 ^ \ 

^' ^ .£1 m + n-(iV + l)/2 ' ^' = ^£, (m + n-(iV + l)/2)2 " ^^^^ 

For the COE we find the form factor 

T-coE^ _ f 2n - n/j()'_„ {S<n<N . 

I '^ I ~ \ 2Ar - n/o^ AT < n , ^^> 

and the variance 

Var(|e'^r)-^co^' = (20) 

r 8n(/i^„ - \f^__^ - i/^+„") - 2n2<_„ < n < iV/2 

8n(/iX. + /\ - \f^,^ - i/^i„") - 2n2^j^_„ + 8iV - 2n iV/2 < n < iV 

I 8n(/o^ + i/o" - i/^+") - 2n2^o^ - 8iV iV < n . 

For the CSE, i.e. the case /3 = 4, Kramers' degeneracy is present. It suffices to take every 
eigenvalue into account once: Assuming the matrix to be of size 2A^ we work with the "trace" 
t^^ = YliLi e~"^* = |tr F" and thus account for every Kramers' doublet once. 
The form factor and its variance are obtained as 

'N/2 




V-n-N/2 0<n<2N 



ie'^p = r:r ^'— ^/^ 2N<n ^^^^ 



Varde'^P) - |tcSE|22 ^ ^22) 

49z:% + t(/:S/. - f~NT) < - < ^ 

-iaZ^T" + !/:"%2 + N-n N<n<2N 

-N 2N <n . 

As the important result we find that fiuctuations are of the same order as the form factor it- 
self. We therefore cannot expect an individual Floquet or random matrix to yield a sequence 
|t„p in accordance with the ensemble mean without any averaging. We shall presently show 
that the form factor is ergodic, i.e. the ensemble mean is approached by taking a local 



time average in Pandey's sense |T0[. Before proceeding to that endeavor we consider the 



fiuctuations in the limit N ^ 00, where they simplify significantly. 



V. ASYMPTOTIC FLUCTUATIONS 

Prior to taking fluctuations to the limit A^ — * oo a normalization is in order. It seems 
quite natural to refer the discrete time n to A^, which is the Heisenberg time for /5 = 1, 2 
and half of the Heisenberg time for (3 = 4:. We write r = n/N ., imagine r kept fixed as 
A^ ^ oo, and eventually allow r to range among the real numbers. The asymptotic form 
factors arising in that limit are of course well known. 

No simplification arises in the unitary case, 



l^-l' -^ K (^\ ] 1^1 < |r| < 1 

-Jf ^^^^^^^ = 11 1 < |r| • (2^) 

The CUE form factor is everywhere continuous but its first derivative jumps at r = and 
at the Heisenberg time, r = ±1. 

In the orthogonal and symplectic cases the sums /^ and g^^ must be turned into integrals 
and obviously yield logarithms. The COE form factor then becomes 



\tn? _ ^ ,. r2|r|-|r|ln(2|r| + l) < |r| < 1 

^^ - KcoEir) = I 2 - |r| ln(|^) 1 < |r| ' (^4) 

We still encounter a jump of Kqq-^{t) at r = but at the Heisenberg time r = ±1 a jump 
arises only for the third derivative. 
The asymptotic CSE form factor. 



l'"l A'cseW^ ?-^'"(l^-™ °1!1^^ (25) 




displays a logarithmic singularity at r = ±1, i.e. at half the Heisenberg time. Jumps arise 
in K'q^-^It) at r = and in -K'cseI''") ^^ ^^e Heisenberg time r = ±2. 

As for the variances Var(|t„p/A^) we find strictly Gaussian behavior in the limit A^ — > oo, 

lim Var(^) = lim I ^^^ 1 = K(rf , (26) 

for all three circular ensembles. A little side remark is indicated, however, for the case of 
the CSE. We here find a funny correction to the Gaussian behavior which according to taste 
one might keep or drop, 

hm Var(M) _ KcsEirf = -5^,1^ , (27) 

where S\r\.i denotes the Kronecker delta. Inasmuch as the r.h.s. for \t\ = 1 (at half the 
Heisenberg time) is of order unity one would want to keep the correction and even feast on 
its non-Gaussian character; on the other hand, if one wants to look at r as a continuous 
variable the exceptional points r = ±1 accomodate removable singularities without weight 
for integrals. Furthermore, the r.h.s. of ( P7| ) is relatively small since i^csE(''") diverges 
logarithmically at r = 1. Figure (|l]) gives an impression for the approach of Var(|t„p/A^) — 



|t„P/A^)^ to the asymptotic — (5|^M7r^/16 as A^ grows. If one wants to study that approach 



analytically, one must isolate in Var(|t„p/A^) — (|t„P/A^)^ the term which does not manifestly 
vanish for A^ — > oo at least as ln(A^)/A^. With the help of (0) and ( ^21) one finds that term 
as 

hin,N) = \-_§li [<l^l (28) 

which can be expressed in terms of the polygamma function ip^^\ the first derivative of the 
digamma function, as 

The polygamma function ip^^^ has the familiar integral representation 

^«(^)= / dt- -e-^K (30) 

JO i — e ^ 

We need the foregoing representation for z = N{1 + 1/(2A^)), z = N{1 — r + 1/(2A^)) and 
z = N{t — 1 + 1/(2A^)). Obviously now, for t 7^ 1 fixed and N -^ 00 we have ip^^^ — ^ in 
all terms involved; at r = 1, however, the limit N —^ 00 yields ip^^\l/2) = 7r^/2 and thus 
h{N, N) = — 7r^/16. The decay of h{TN, N) to zero for A^ — > 00 at fixed r 7^ 1 is found from 
the asymptotic large- 2; behavior of the polygamma function 



1 „, 1 

2- 



z z- 

as 



^«(^) = ^ + 0(-^) (31) 




h{rN^N) = \ BA^V-3^ °<^<^ _ (32) 



VI. ERGODICITY OF THE FORM FACTOR 

Fluctuations of the form factor tend to be suppressed by a local time average, 

^ n+An/2 

(|tnH = ^ E l^n'P. (33) 

n'=n— An/2 

As one increases the time window An one expects the time average to become equivalent to 
the ensemble average. Inasmuch as the systematic time dependence should not be washed 
out by the time average, one must first let A^ — ^ 00 and subsequently An -^ 00. 

To uncover the expected ergodicity we consider the form factor as normalized to its 
ensemble mean, |tnP/|^nP- We realize that (|tnP/|^nP) = 1 holds trivially and propose to 
show that the ensemble variance of the temporal mean vanishes, Var((|t„p/|t„p)) -^ 0, as 
An — »• 00. The variance in question is defined as 



Var (|t„|V|tnP) 



1 \2 n+An/2 / ,, pn p ^ 

^^^ n',n"=n-An./2 Vl^n'P |^n"P y 



(34) 



Its evaluation requires knowledge of the cross-correlator |tn'n^n"P for n' ^ n" . We may 
determine the cross-correlator from the characteristic function 



P{kn', K") = exp ( -- ^ ((A;„,e-'«V. + fc„„e-i"'V.) + c.c.) 



which generalizes the single-trace one (^) and yields 



itn'nt„»p = (2ir 



dkn/dkl,dkn"dk*^/, 



P{kn',kn") 



(35) 



(36) 



Once more employing the technique used to establish the single-trace characteristic function 
one finds for the CUE 



P'^{kn',kn") =detM 



(37) 



M, 



flU 



'" ^pi{M--)</'exp ( -^(A;,,e-'"'* + A;„»e-^""<^ + c.c.; 



27r 



and from here (n' > n") 



\f /Pit ,,|2 
I ^n I I *^n I 

17 [2 17 [2 



N-n'-n" 

■n'r)" 



Af-n'+n" 







Nn" 



n' < N, n' + n" > N 
N < n', n' -n" < N 
otherwise . 



(38) 



Upon inserting these cross-correlations into ( ^41) the desired variance is obtained. An upper 
bound for this variance is found by neglecting all negative contributions to the sum, i.e. all 
terms with n' ^ n" . The estimate then reads 



Var((|t„|VK:F))<-^^^°°0. 



An 
The asserted ergodicity of the form factor is thus proven. 



(39) 



VII. ASYMPTOTIC INDEPENDENCE OF THE TRACES 



With the cross-correlations ( PBD at hand we can now show that the traces t^'j^n" with 
n' 7^ n" are asymptotically uncorrelated in the limit N ^ oo. Again we assume n' > n" and 
first consider the case n' < N, n' + n" > N m (RSI) , where we can estimate 



N -n' - n" 



< 



n 



n'n" n' N 



(40) 



In this situation n' cannot be small compared to A^ since n' > n" and n' + n" > N. In the 
case N < n'n' — n" < N the estimate reads 



N -n' + n" 



Nn" 



< 



n 



Nn" N 



(41) 



We see that the cross-correlations (|38D vanish no more slowly than 1/A^ for A^ — > cxo. 



VIII. THE FORM FACTOR OF THE KICKED TOP 

The Floquet operator of the kicked top with a chaotic classical limit has proven faithful 
to RMT predictions in many aspects [|l|. Among these, the level-spacing distribution and 
the integrated two-point correlator are noteworthy. We here propose to reveal faithfulness 
of the top to RMT for the form factor and its fluctuations. 

The Floquet operator is built of angular momentum operators J^, Jy, Jz and can be 
understood as a succession of linear rotations and torsions. For integer angular momentum 
values j the operator to be considered is of the form 

F = Q-^2j+i-^^-^I^^J^Q-^l^yJyQ-^2]+i'^S-^l^^'^^ (42) 

and depending on the rotation angles /5j and torsion constants Tj belongs to the orthogonal 
{tz = 10, Pz = l,Py = l,r^ = Px = 0) or unitary (r^ = 10, /3^ = l,/?^ = l,r^. = 4,/3^ = 1.1) 
universality class. For semi-integer j the Floquet operator 

ri = l,r2 = 4,r3 = 2.1,r4 = 10 (43) 

lacks all geometric symmetries and thus belongs to the symplectic universality class [|1|. 
For numerical work F is represented as a matrix of dimension A^ = 2j ' + 1 in the basis of 
eigenvectors of J^, Jz\jfT^) = m\jm), with fixed j and —j<m< +j. 

Figure § shows the form factors for the three cases chosen and reveals wild fluctuations 
about the ensemble means which are displayed as the solid lines. The dashed lines together 
with the abscissa border stripes of one and two asymptotic linear variances around the means. 
Gaussian distributions for the t„ would predict the fractions (1 — l/e)/w(e) ^ 0.383 and 
(1 — 1/e^) ~ 0.865 of all traces to lie within the one-variance and the two-variances stripes, 
respectively. On counting we find the fractions to be 0.390 and 0.875 for the orthogonal top, 
0.397 and 0.870 for the unitary one, and finally 0.386 and 0.858 for the symplectic top. The 
Gaussian expectation is borne out well by these numbers and further substantiated by the 



histograms in Fig. ^ for the normalized moduli r„ = \tn\ / v Itn]"^ ■ Gaussian statistics would 

2 

yield the distribution 2re~'^ which is displayed as the smooth line and well approximated 
by the histograms. 

Finally, we illustrate ergodicity of the form factor of the top (unitary variant) in Fig. ^. 
The left part shows the form factor fluctuations after a local time average over the window 
An = 20. Comparison with Fig. |^ reveals the smoothing effect of the local time average. We 
have also evaluated the variance of the "band" of points It^P/ltnP as a function of the time 
window An. The result is represented by the dots in the right part of Fig. |^, together with 
the 1/An bound predicted by RMT. Again the faithfulness of the top to RMT is impressive. 
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FIGURES 



FIG. 1. The approach of the symplectic ensemble's Var(|t„|^/A^) — (|t„p/A^)^ to the asymptotic 
— (5|^M7r^/16 with growing N at half of the Heisenberg time. 

FIG. 2. Form factors of the orthogonal, unitary and symplectic Floquet operators with param- 
eters as given in the text. Also shown are the respective ensemble means and stripes of one and 
two linear variances width around the mean. 



FIG. 3. The distributions of the Floqet operators' normalized traces r„ = |in|/\/ l^np displayed 

2 

by the histograms show good agreement with the distribution 2re~'^ (smooth line) expected for 
Gaussian statistics. 

FIG. 4. Left: Form factor of the unitary top smoothened by a local time average over the 
window An = 20. Right: Width of the fluctuating band for |tnP/|inP s^s a function of the time 
window An (dots) together with the 1/An bound (line) predicted by RMT. 
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